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Abstract. 

We  consider  interactions  between  varicose  and  sinuous  oblique  disturbances 
in  the  Bickley  jet,  using  both  nonlinear  stability  theory  (in  its  nonlinear 
critical  layer  form)  and  direct  numerical  simulation  using  a spectral  method. 
Nonlinear  stability  theory  indicates  that  a (nonlinear)  interaction  between 
the  modes  should  occur,  and  our  simulations  would  seem  to  support  this. 


1.  Introduction 

This  study  is  concerned  with  three-dimensional  transition  in  plane  wakes 
and  jets.  Nonlinear  stability  theory  [9,  10]  has  predicted  that  modal  inter- 
actions play  an  important  role  in  that  transition,  and  in  this  paper,  we  will 
first  present  an  overview  of  that  theory  and  then  present  simulations  that 
support  the  theory. 

For  the  plane  wakes  and  jets  considered  here,  it  is  well  known  that 
there  may  be  two  different  types  of  neutral  modes  with  critical  layers  cen- 
tred on  the  inflection  points,  viz.  the  sinuous  and  varicose  modes.  The 
plane  (Bickley)  jet,  which  has  a sech2y  velocity  profile  has  been  used  by 
numerous  authors  to  provide  a good  approximation  to  such  a wake  behind 
a bluff  body,  is  somewhat  special  in  that  the  varicose  and  sinuous  modes 
have  neutral  wavenumbers  of  1 and  2 respectively,  so  that  the  former  is 
the  subharmonic  of  the  latter.  Several  studies  have  explored  the  possibility 
of  an  interaction  between  these  two  modes,  and  this  sort  of  interaction  is 
the  focus  of  the  present  study.  For  two-dimensional  disturbances,  Kelly  [6] 
used  Stuart- Watson  type  nonlinear  stability  theory  to  investigate  interac- 
tions of  this  type;  however,  he  found  that  there  was  no  modal  interaction  of 
the  type  assumed.  Later,  Leib  and  Goldstein[8]  re-examined  the  problem  for 
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purely  two-dimensional  disturbances  using  a nonlinear-nonequilibrium  crit- 
ical layer,  and  they  found  that  there  was  indeed  an  interaction  between  the 
modes.  Mallier[9,  10]  studied  three-dimensional  disturbances,  since  Gold- 
stein]^, 2]  had  earlier  shown  that  a pair  of  oblique  waves  superimposed 
on  a shear  layer  could  interact  nonlinearly  to  give  rise  to  extremely  rapid 
growth,  and  also[4]  that  when  the  oblique  waves  were  inclined  at  ±60°, 
an  additional  interaction  could  take  place  between  the  oblique  waves  and 
a plane  wave,  so  that  the  growth  was  faster  still;  this  last  mechanism  is 
known  as  a “resonant  triad”.  Mallier[9]  explored  the  possibility  of  an  in- 
teraction between  a pair  of  resonant  triads  in  the  Bickley  jet,  with  one 
triad  consisting  of  a plane  sinuous  mode  together  with  a pair  of  oblique 
sinuous  modes  inclined  at  ±60°  and  the  other  triad  consisting  of  a plane 
varicose  mode  together  with  a pair  of  oblique  varicose  modes  also  inclined 
at  ±60°;  the  motivation  behind  this  was  to  see  if  the  triads  could  interact 
so  that  the  growth  was  more  rapid  than  if  only  a single  triad  were  present. 
Mallier  found  that  interactions  could  occur,  and  his  study  essentially  cov- 
ered three  stages:  a linear  stage  when  the  amplitudes  of  the  disturbances 
were  very  small,  the  “parametric  resonance”  stagefll],  and  the  so-called 
“fully-coupled”  stage[4,  21,  12].  The  amplitude  equations  presented  were  of 
course  for  the  third  (fully-coupled)  stage,  but  the  two  earlier  stages  could 
be  recovered  from  these  equations  by  rescaling  the  amplitudes,  as  discussed 
in[4].  The  study  of  the  fully-coupled  stage  was  a little  restrictive  in  that 
it  was  necessary  to  assume  that,  in  that  stage,  the  varicose  oblique  modes 
were  larger  any  of  the  other  waves  present  which  is  at  odds  with  the  lin- 
ear theory  which  says  that  the  linear  growth  rates  of  the  sinuous  modes 
are  larger  than  those  of  the  varicose  modes.  This  was  because  in  the  earlier 
parametric  resonance  stage,  when  it  was  assumed  that  all  of  the  waves  were 
of  the  same  order  of  magnitude,  it  had  been  found  that  the  varicose  oblique 
waves  underwent  very  rapid  growth  while  the  plane  waves  and  the  sinuous 
oblique  waves  continued  to  grow  exponentially  in  a linear  fashion.  In  ad- 
dition, the  coupling  in  the  equations  in  [9]  was  a little  unusual  in  that  the 
sinuous  triad  did  not  affect  the  varicose  triad,  and  therefore  equations  for 
the  varicose  triad  were  simply  those  for  a single  resonant  triad[4,  21,  12]. 
However,  the  sinuous  triad  was  strongly  affected  by  presence  of  the  varicose 
triad,  and  furthermore,  if  the  varicose  triad  was  absent  the  nonlinear  terms 
in  the  equations  for  the  sinuous  triad  vanished,  leaving  only  linear  equations 
for  those  modes.  Mallier [10]  later  studied  the  case  of  two  pairs  of  oblique 
waves  superimposed  on  the  Bickley  jet  at  the  same  angle,  ±9:  one  pair  was 
varicose,  the  other  sinuous.  Once  again,  an  interaction  was  found  to  occur 
between  the  modes.  In  both  studies,  [9,  10]  the  end  result  was  a set  of 
highly  nonlinear  coupled  (Hickernell-type[5])  integro-differential  evolution 
equations,  the  solutions  to  which  had  a finite-time  singularity.  These  equa- 
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tions  involved  a nonlinear  kernel,  with  the  nonlinearity  being  cubic  in[10] 
and  quartic  in[9] . Both  studies[9,  10]  used  nonlinear  critical  layer  theory 
and  followed  the  approach  taken  earlier  by  Goldstein[3,  4]  for  a flow  with 
a single  critical  layer  and  a single  unstable  mode. 

The  reason  these  interactions  between  the  varicose  and  sinuous  modes 
are  considered  important  is  that  it  is  possible  they  can  cause  extremely 
rapid  nonlinear  growth.  Unfortunately,  at  the  time  the  studies  mentioned 
above  were  performed,  there  was  little  if  any  experimental  or  numerical 
evidence  to  corroborate  our  analysis,  and  the  situation  remains  the  same 
today.  Some  experiments  have  hinted  at  interactions,  but  have  not  explored 
it  further.  Wygnanski  et  al.  [23]  conducted  careful  experiments  on  small 
deficit  (turbulent)  wakes  and  found  that  the  development  of  some  aspects 
of  the  flow  was  dependent  on  initial  conditions,  which  they  attributed  to 
interactions  between  the  varicose  and  sinuous  modes,  and  other  experiments 
(e.g.[13,  14])  have  also  suggested  that  these  interactions  may  take  place. 
We  should  also  mention  that  very  rapid  amplification  of  three-dimensional 
disturbances  has  indeed  been  observed  in  plane  wakes  in  both  experiments 
(e.g.  [1,  19,  20])  and  numerical  simulations[17,  18],  but  it  is  unclear  (at 
least  to  the  present  author)  how  much  of  that  growth  is  attributable  to  the 
Goldstein  mechanisms  for  three-dimensional  instability  and  how  much  is 
due  to  an  interaction  between  the  varicose  and  sinuous  modes. 

In  an  attempt  to  remedy  what  we  perceive  as  a lack  of  numerical  veri- 
fication of  the  theoretical  results,  we  present  here  some  preliminary  results 
from  direct  numerical  simulations  of  the  Bickley  jet  and  compare  those  re- 
sults to  the  predictions  of  our  earlier  asymptotic  analysis.  The  outline  of 
the  rest  of  the  paper  is  as  follows.  In  the  next  section,  we  review  the  theory 
and  resulting  amplitude  equations.  After  that,  we  give  the  details  of  our 
numerical  method  and  present  the  results  of  our  simulations.  Finally,  we 
make  some  concluding  remarks. 

2.  Review  of  Theory 

In  [9,  10],  we  considered  the  stability  of  the  Bickley  jet  u0  = (sech2y,0, 0^ 
to  three-dimensional  disturbances,  either  in  the  form  of  two  pair  of  oblique 
waves, 

2 eAnuueiacos0(x-ct)  cos  [asinflz]  + c.c. 

+ 2sA22u22e2tacose^  cos  [2a sin 0*]  + c.c.  (1) 

or  in  the  form  of  a pair  of  resonant  triads,  consisting  of  a two  pair  of  oblique 
waves  at  ±60  together  with  two  plane  waves, 

2eA\\ullelct^x~ct^‘2,  cos  a\/3z/2j  + c.c. 
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+ 2eA22U22e2ia(x-ctV2  cos  [aV3z]  + c.c. 

+ e4/3  [a2 ow20eiQ(l“ci)  + A4 oM4oe2ia(l~ct)  + c.c.]  . (2) 

Prom  the  linear  theory,  we  know  that  the  Bickley  jet  has  two  inflection 
points  at  y = ±arccosh  i/3/2  where  no  = 2/3,  and  also  two  neutral  modes 
with  c = 2/3  which  have  critical  layers  centered  on  the  inflection  points:  a 
varicose  mode  with  a = 1 and  v = sech  y tanh  y [15]  and  a sinuous  mode 
with  a = 2 and  {>  = sech2y  [16].  In  the  above  disturbances,  we  have  either 
two  amplitudes  (for  the  pairs  of  oblique  waves)  or  four  amplitudes  (for  the 
resonant  triad)  and  the  objective  of  [9,  10]  was  to  derive  amplitude  equa- 
tions for  those  amplitudes,  and  examine  how  they  affected  each  other.  In  the 
analysis,  it  was  assumed  that  the  modes  were  periodic  in  time  and  spatially 
growing,  with  the  wavenumber  a taking  the  neutral  value  of  1 and  phase 
velocity  being  perturbed  slightly  from  neutral  c = 2/3— nc\.  This  perturba- 
tion from  neutral  led  to  a slow  time  scale  T = /it  on  which  the  amplitudes 
evolved.  For  extremely  small  disturbances,  it  was  found  that  the  growth  was 
linear,  with  Au(T)  = A$eanX  and  similar  expressions  for  the  remaining 
modes,  but  for  larger  disturbances  the  evolution  was  nonlinear;  the  evo- 
lution first  became  nonlinear  when  e = /A  For  pairs  of  oblique  waves,  it 
was  found  that  the  amplitude  equations  were  (Hickernell-type)  integrod- 
ifferential  equations  with  a cubic  nonlinearity.  For  the  resonant  triad,  the 
equations  were  of  a similar  form  but  with  a quartic  rather  than  a cubic  non- 
linearity. In  [10],  we  solved  these  equations  numerically,  using  Goldstein’s 
numerical  scheme  [3],  and  it  was  found  that,  as  with  similar  problems,  the 
evolution  of  the  disturbances  went  through  3 stages:  initially,  the  distur- 
bances grew  linearly,  until  a second  finite-amplitude  nonlinear  stage  was 
reached,  and  eventually,  the  oblique  waves  experienced  explosive  growth. 
Goldstein  [3]  showed  that  his  equations  had  a singularity  after  a finite  time 
Ts  (or  at  a finite  distance  downstream),  and  was  able  to  fit  a structure  to 
it  A ~ ao(Ts  — T)-3"^.  This  same  structure  applies  to  both  modes  for 
the  Bickley  jet.  Although  our  study  did  not  include  viscous  effects,  other 
studies  for  related  problems  [7,  22]  have  shown  that  weak  viscosity  can 
delay  the  onset  of  this  finite  time  singularity  but  not  eliminate  it,  so  that 
our  results  are  still  meaningful  at  high  Reynolds  numbers.  The  origins  of 
this  finite-time  singularity  are  still  not  entirely  clear,  but  it  appears  to  be 
connected  to  the  breakdown  of  the  theory  and  the  onset  of  a new,  still  more 
nonlinear  stage  governed  by  the  full  Euler  equations. 

3.  Numerical  Simulations 

We  now  turn  to  our  numerical  simulations,  in  which  we  have  employed  a 
standard  spectral  (Fourier)  method.  In  our  computations,  we  decomposed 
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Figure  2.  3D  Run2 

the  velocity  into  a mean  flow  and  a perturbation,  u ~ (uo(y),0, 0)  + u,  and 
then  assumed  that  the  mean  flow  is  independent  of  time  and  the  viscosity 
acts  only  on  the  perturbation;  in  reality,  this  would  require  the  presence 
of  a body  force  to  counteract  the  effect  of  viscosity  on  the  mean  flow.  The 
perturbation  therefore  obeys 

du  1 

— = —u  • Vu  — Vp  + — V u 
ot  Re 

V*u  = 0.  (3) 

The  approach  we  took  was  to  assume  that  the  flow  was  periodic  in  both  x 
(streamwise)  and  2 (spanwise),  and  use  complex  Fourier  series  in  those  di- 
rections. We  truncated  the  y-direction,  so  that  our  domain  was  — Y <y<Y 
rather  than  — oo  < y < oo  and  used  sines  and  cosines  in  that  direction.  We 
calculated  the  nonlinear  terms  in  physical  space  and  derivatives  in  Fourier 
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Figure  4-  3D  Run4 


space,  and  used  a fast  fourier  transform  (FFT)  to  switch  between  real  and 
physical  space.  In  our  computations,  it  was  necessary  to  calculate  the  pres- 
sure, which  we  were  able  to  do  by  using  incompressibility  and  then  inverting 
a Laplacian, 

p = — (v2)  _1  V • jVX;  (4) 

this  is  fairly  easy  to  do  with  a Fourier  method.  For  the  time-stepping,  we 
used  an  explicit  Adams-Bashforth  scheme  for  the  nonlinear  terms, 

ui  = u0  + — (3 F0  — F_i) , (5) 

while  for  the  viscous  terms,  a semi-implicit  Adams-Moulton  scheme  was 
used, 


ui  = uq  -(-  — ( Fi  + Fq)  . 

In  our  simulations,  we  calculated  the  energy  in  each  mode, 


I umn\2  dzdxdy, 


(6) 

(7) 


and  compared  this  to  the  theoretical  amplitude  of  the  disturbance.  We 
should  mention  however  that  strictly  speaking,  therefore,  we  are  not  com- 
paring like-with-like,  since  the  two  definitions  differ  slightly.  Similarly,  our 
simulations  include  weak  viscosity  (a  Reynolds  number  of  1500  was  used) 
which  is  necessary  for  numerical  stability,  while  our  theory  was  inviscid. 
However,  as  we  discussed  earlier,  studies  by  other  authors  have  indicated 
that  our  theory  is  applicable  to  high  Reynolds  number  flows. 


3.1.  SAMPLE  3D  RUNS 

In  Fig.s  1-4,  we  present  several  different  runs  for  the  resonant  triad,  with 
6 — 7t/3.  Each  of  the  runs  shown  eventually  blows  up  as  we  lose  spectral 
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decay  and  the  higher  harmonics  become  too  large.  Paradoxically,  this  is 
similar  to  the  reason  for  the  finite-time  singularity  in  Goldstein’s  amplitude 
equation.  Since  the  location  of  the  blow-up  is  different  for  each  of  the  runs 
presented,  this  by  itself  would  indicate  that  an  interaction  between  the 
modes  is  occuring.  Our  simulations  indicate  that  we  are  able  to  capture  the 
nonlinear  stage  that  appeared  in  the  solution  to  Goldstein’s  equation,  and 
for  the  three-dimensional  case  it  would  appear  that  the  modes  do  indeed 
exert  a significant  influence  on  each  other;  an  example  of  this  can  be  seen 
in  the  behavior  of  A\\  in  Fig.s  1 and  2:  the  initial  conditions  for  An  were 
the  same  in  both  runs,  and  the  difference  in  the  behaviors  of  this  mode  in 
the  two  runs  is  due  entirely  to  nonlinear  interactions  between  the  modes. 
The  initial  conditions  for  the  runs  shown  are  as  follows:  in  Run  1,  A20  and 
A40  were  initially  zero;  in  Run  2,  A20,  A40  and  A 22  were  all  originally  zero; 
in  Run  3,  A20,  A40  and  An  were  all  originally  zero;  in  Run  4,  A40  was 
all  originally  zero.  One  point  to  notice  about  Run  3 is  that  since  An  and 
.A20  were  both  initially  zero,  they  remain  zero:  those  two  modes  cannot  be 
generated  by  the  interaction  of  the  other  two  modes.  It  is  interesting  to  note 
in  these  figures  that  the  behaviour  of  one  mode  depends  upon  which  other 
modes  is  present,  which  confirms  that  there  is  indeed  a strong  nonlinear 
interaction  between  the  modes  as  suggested  by  [9,  10]. 

4.  Conclusions 

In  the  preceding  sections,  we  have  outlined  the  theory  presented  in  [9,  10] 
for  nonlinear  interactions  in  the  (plane)  Bickley  jet  between  the  varicose 
and  sinuous  modes,  and  then  presented  some  of  the  three-dimensional  re- 
sults obtained  using  DNS.  For  both  the  three-dimensional  case  and  also 
the  two-dimensional  case  [8],  the  nonlinear  theory  suggests  that  there  is 
an  interaction  between  the  sinuous  modes  and  the  varicose  modes.  The 
three-dimensional  simulations  presented  here  would  appear  to  confirm  that 
because  the  behaviour  of  a mode  clearly  differs  depending  upon  which  other 
modes  are  present.  This  behaviour  is  quite  strong  in  the  three-dimensional 
case  and  two-dimensional  simulations  not  presented  here  indicate  that  the 
two-dimensional  modes  also  interact  but  that  the  interaction  in  that  case 
is  much  weaker.  We  are  currently  performing  more  simulations,  although 
this  is  a slow  process  (the  turnaround  time  for  3D  runs  on  the  UWO  Cray 
is  about  2 months),  and  we  hope  to  present  more  complete  results  at  some 
point  in  the  future. 
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